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>and  an  increase  in  the  field  at  the  wall.  Reflection  here  means 
re-emission,  i.e.,  a  return  as  half-Maxwellian.  < 

Single  ended  plasma  device,  general  behavior  dc  or  ac^-This  lausanne 
^  --IOPP  paper  describes  simulations  with  various  dc  biases,  showing  both 
stable  states  (floating,  negative  bias)  and  unstable  (moderate  positive 
bias)  with  periodic  current  disruptions  (relaxation  oscillations) 
for  both  electron  and  ion  rich  sources.  '  ^ '' 

Jh  Single  ended  plasma  device,  unstable  states*-- This  Innsbruck  Double 
„  _Layer  Conferences’ paper  is  more  detail  of  the: instabil ities  in  0,  ''  . 

showing  J(t)  and  phase  space  for  the  ion-rich  instability,  (t  ‘ 

£•.  Corrections  to  time-independent  Q-machine  equil ibria^-In  the  ion-rich 
source,  negatively  biased  collector  simulations,  the  classical  analysis 
has  an  electron  hole  at  the  source.  However,  this  hole  is  filled 
in  the  simulations,  even  when  the  number  of  particles  is  made  very 
large  (reducing  fluctuations).  The  corresponding  filled-hole  analysis 
fits  the  measured  potentials.  •'  - 

F.  Multifluid  derivation  of  the  Alfven  ion-cyclotron  linear  dispersion 
relation«fThis  is  an  alternative  derivation  which  allows  direction 
connection^ between  features  of  the  wave  responsible  for  the  currents 
and  terms  in  the  dispersion  relation. 

G.  Potential  barrier  between  hot  and  cool  plasmas.  These  simulations  are 
aimed  at  explaining  the  experimental  observation  of  a  potential  dip  be¬ 
tween  hot  and  cool  plasmas.  So  far  no  spatial  interface  has  been 
observed. 

Code  Development 

.  A.  P0M1 ,  Particle  simulation  of  non-periodic  (bounded)  systems.  The  code 
is  out,  with  a  User's  Manual',"  in  heavy  use  here  (see  A,B,C,D,E,G 
above)  and  available  to  others. 
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SECTION  I:  PLASMA  THEORY  AND  SIMULATION 


A.  ELECTRON  BERNSTEIN  WAVE  INVESTIGATIONS 

Wm.  S.  Lawson  and  J.  P.  Lynov  (visitor,  Risa,  Denmark) 


A  modification  of  the  code  PDWl  is  being  used  to  investigate  the  properties 
of  oblique  Bernstien  waves  at  finitesimal  amplitudes.  So  far,  the  code  has  been 
shown  to  work,  and  reproduce  Bernstein  waves  with  the  correct  damping  rate 
(for  oblique  angles).  A  significant  wave  number  shift  has  been  noted,  however, 
for  relatively  small  excitations. 

The  numerical  experiment  is  modeled  on  the  Double  Plasma  experiment  of 
Armstrong,  Rasmussen.  Stenzel,  and  Trulsen  [l].  The  exciter  is  composed  of 
three  grids,  the  outer  two  being  tied  together,  and  being  voltage  driven  with 
respect  to  the  center  grid.  This  experimental  model  may  be  responsible  for 
some  peculiar  effects  observed  at  Low  amplitudes,  such  as  the  apparent 
reflection  of  the  driven  wave. 

Other  effects  observed,  but  not  explained  are  the  apparent  excitation  of 
waves  at  higher  frequencies  than  the  driving  frequency,  and  reduced  damping 
at  high  amplitude.  Work  is  continuing  in  an  attempt  to  explain  these  results. 

An  ongoing  problem  is  the  relatively  high  noise  level  of  the  simulations. 
This  noise  seems  to  be  associated  with  the  high  dimensionality  of  the  simulation 
(1-D.  3-V),  and  has  dictated  large  numbers  of  particles  and  time  averaging. 

[l]  Armstrong  et.  al.,  Physics  Letters  85A(i981)291-.284. 


B.  SIMULATION  OF  THE  PLASMA-SHEATH  REGION 

INCLUDING  ION  REFLECTION 

Lou  Ann  Schwager  (Prof.  C.  K.  BirdsaTI,  Dr.  I.  Roth) 


A  low  temperature  plasma  interacts  with  a  collector  plate  primarily  through  charged  particle  reflec¬ 
tion, emission,  or  absorption.  These  phenomena  will  affect  the  electric  potential  profile  and  hence  the  trans¬ 
port  to  the  plate.  Present  simulation  studies  indicate  that  a  reflected  ion  current  increases  the  plasma 
potential  above  that  simulated  with  a  purely  absorbing  plate. 

The  computer  code  PDWl  models  this  behavior.  The  setup  of  this  simulation  is  identical  to  that 
described  in  the  previous  report  (QPR  III, IV  1983).  This  earlier  study  modelled  plasma  flowing  into  a  purely 
absorbing  plate.  The  current  work  differs  from  the  previous  only  in  the  inclusion  of  ion  reflection  at  the 
plate.  The  system  is  shown  in  Figure  1.  As  in  the  earlier  run,  ions  and  electrons  returning  to  the  source  are 
thermalized  and  re-injected  with  the  initial  flux.  We  have  chosen  a  maw  ratio  of  40  to  conserve  computation 
time  and  an  ion  reflection  coefficient  of  70%.  This  coefficient  represents  those  parameters  which  are  typical 
of  the  plasma  region  at  the  collector  plate  in  the  end  cell  of  a  tandem  mirror  reactor  (and  possibly  in  the 
poloidal  divertor  of  a  tokamak). 

For  a  hydrogen  plasma  with  temperatures  of  30-100  eV  incident  on  a  plate  of  molybdenum,  the  atomic 
data  in  Table  1  pertains. 

Table  1.  Atomic  Data  for  30-100  eV  H  on  Mo 


Mechanism 

Reflected  Ion  * 
per  Ion 

Secondary  Electron 
per  Electron 
Reflected  Electron 
per  Electron 
Secondary  Electron 
per  Ion 

Sputtered  Mo  atom 
per  Ion  (1  keV) 


Interaction 
Coefficient 
.6-.8  1 

.4-.8  2 

.05-.  10  2 

0 2 

2xl0“3 2 


*  It  was  later  discovered  that  only  a  small  fraction  of  the  reflected  particles  are  ionized.  This  will  be  accounted 
for  in  later  reports.  _ 


Obviously  the  plate  does  not  behave  as  a  purely  absorbing  boundary.  We  model  ion  reflection  here 
because  it  dominates  the  other  mechanisms.  In  future  work  we  intend  to  include  others  from  Table  1. 

Unfortunately,  little  information  is  available  on  the  reflection  of  low  energy  light  ions.  Oen  and  Robinson 1 
used  the  simulation  program  MARLOWE  to  derive  universal  curves  for  the  reflection  coefficients  of  various 
light  ions  independent  of  target  material.  They  plot  the  reflection  coefficient  as  a  function  of  the  Thomas- 
Fermi  reduced  energy,  e.  The  formula  they  use  is 

0.03  EM  3 

€  “  zZ{m  +  M )(z3'3  +  Z2/3)i/2 

where  m  and  M  are  the  projectile  and  target  atomic  weights,  z  and  Z  are  the  projectile  and  target  atomic 
numbers,  and  E  is  the  projectile  energy  in  eV.  Hence,  for  energies  of  30-100  eV,  protons  hitting  a  Mo  target 
have  e  ranging  from  0.02  to  0.006.  With  this  range  for  £,  Figure  6  of  Reference  2  specifies  that  the  ion 
reflection  coefficient  lies  between  60%  and  80%. 

Electric  potential  histories  at  x  =  L/2  for  reflection  coefficients  of  0%  and  70%  are  shown  in  Figures  2A 
and  2B.  Electrons  reaching  the  wall  before  ions  cause  the  early  potential  rise  in  both  cases.  The  interaction  of 
the  opposing  ion  flows  prior  to  equilibration,  seen  in  Figure  3A,  generates  the  perturbed  potential  profile  in 
Figure  3B.  This  appears  as  a  mild  transient  in  the  potential  in  Figure  2B.  After  equilibration,  the  opposing 
ion  flows  lead  to  the  mildly  perturbed  potentials,  as  shown  in  Figure  4A.  This  results  in  the  nearly  flat 
potential  profile  of  Figure  4B.  Note  also  in  Figure  3B  that  the  70%  case  generates  a  steeper  potential  profile 
in  the  sheath  region  (where  Ewau  ~  —19)  than  that  for  the  0%  case  in  Figure  4C  (where  Ewau  =  —5.8). 
The  steep  potential  rise  away  from  the  wall  repels  low  energy  reflected  ions  at  x  =  0.  The  ion  excess  in  the 
sheath  then  generates  an  even  stronger  electric  field. 

Our  next  steps  are  to  include  the  study  of  a  reflected  flux,  which  depends  on  the  incident  flux  and 
energy,  and  of  the  effects  of  secondary  electron  emission. 

References 

1.  O.S.Oen  and  M.T.Robinson.  "Computer  Studies  of  the  Reflection  of  Light  Ions  from  Solids”,  Nuclear 
Instruments  and  Methods  132,  pp.647-653  (1976). 

2.  C.F. Barnett  et  aL,  "Atomic  Data  for  Controlled  Fusion  Research”,  Oak  Ridge  National  Laboratory. 
Oak  Ridge,  Tennessee,  ORNL-5207,  pp.D.1.8-D.5.13  (1977). 

3.  T.J. Dolan,  Fusion  Research,  Vol  HI,  Pergamon  Press,  New  York,  1980,  p.697. 


-3- 


CO 


O 

Z 

O 

LkJ 

CO  -J 

Z  UJ 

s  « 

co  to 

O  2  Z 

“  o  o 


UJ 

o 

UJ 

H 

O 

UJ 

-a 

Z 


Q 
UJ 
X 
D 
— J 
U. 
UJ 

z 


Q 

UJ 

I- 

o 

UJ 


^  U 
2° 

_iO 
0-  CO 


CO 

CO 

Id 

z 

z 

o 

CO 

o 

o 

h- 

UJ 

z 

'  < 

(D 

z 

o 

to 

as 

o 

z 

& 

UJ 

_J 

a 

UJ 

h 

U 

UJ 

a 

UJ 

CD 

z 

-J 

CL 

< 

UJ 

-J 

o 

u. 

CO 

UJ 

GO 

z 

< 

-4- 
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E=0 


SYSTEM  LENGTH 
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Figure  4.  Profiles  in  (a)  ion  pnase  space 
for  707.  reflection,  (b)  electric  potential 
for  70%  reflection,  and  (.c)  electric  poten¬ 
tial  for  0%  reflection.  All  at  the  last  time 
step  (at  equilibrium) . 
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Figure  1.,  parameter  apace  for  the  aingle-ended  collisionless  Q-machine, 
a  ▼».  rjep;  where  a  =  n+  /n+  the  neutralization  parameter,  and  rjcp  = 
e^c pJkT  the  normalized  wall  potential  or  bias. 


/£  =  exp(-i/J)  •  H(u,rj(X)) 


where  n  =  m*/Mi  (and  we  have  normalized  to  the  density  of  the  source 
electrons, 

The  two  parameters  a  and  rjep  (the  normalized  cold  plate  bias)  define 
the  parameter  space  for  the  single-ended  Q-machine,  shown  in  Figure  1. 
We  are  going  to  consider  only  the  single  maximum  potential  profile,  which 
occurs  in  the  upper  left  region  of  parameter  space,  negative  bias  with  ion 
rich  emission(Le.  a  >  l).  Figure  2  shows  the  assumed  potential  profile 
for  the  analysis.  We  now  model  the  distribution  functions  as  functions  of 
potential.  For  the  electrons,  for  rj  >  0  (x  <  xa)  we  have 

/«+(«'.»?(x))  =exp(-(i/3  -»Kx)))  '  Hil I'-n/SEx)) 

/.~(*/.»?(x))  =  oq>(-(x/2-i7(x)))  •  B(v+y/ri{x))  •  H{y  -  yjnix)  -  rjcp) 


where  we  have  denoted  the  positive  and  negative  part  of  the  distribution 
function  by  +  and  -  respectively.  The  low  velocity  cutoff  for  /+  is  due  to 
the  fact  that  the  electrons  are  initially  accelerated  by  the  potential  bump  as 
they  enter  the  system  from  the  left.  The  same  is  true  of  the  term  in  f~  for 
the  returning,  reflected  electrons  (i.e.  in  the  potential  bump  the  minimum 
energy  a  particle  should  have  is  1/2 m«s  =  e<f>( x)).  This  leads  to  a  *hole’  in 
electron  phase  space,  which  can  be  seen  in  Figure  2  where  we  have  drawn 
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distribution  characterised  by  the  temperature  of  the  hot  plate.  The  elec¬ 
trons  follow  the  Richardson-Dushman  equation  for  thermionic  emission, 
this  leads  to  equal  electron  and  ion  temperatures  (r,  =  Te  =  T^otpiatc).  We 
then  model  the  distribution  for  each  species  at  the  source  by, 

0-.+  -.2 


where  vt,  =  (2 kT/m,)1/2t  and  k  is  Boltzmann’s  constant,  and  a  denotes  the 
species  (i.e.«  =  e,»),  and  H(v)  is  the  Heaviside  step  function.  For  conve¬ 
nience  we  now  introduce  the  following  normalisation, 


x 


_  v 
~  VU 


( length ) 
(velocity) 

( potential ) 


where  XDe  =  (KT/4icneoe2)1/2  is  the  Debye  length  associated  with  the  elec¬ 
tron  density  of  the  source,  and  define  the  neutralization  parameter  V  as, 


a  =  2*o 

nto 


the  ratio  of  the  ion  to  the  electron  density  leaving  the  source. 

We  may  now  write  the  distribution  functions  for  the  source  ions  and 
electrons  in  normalized  form  as 
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E .  Corrections  to  Time  Independent  Q-Machine  Equilibria 
Perry  C.  Gray  and  Thomas  L.  Crystal 

The  ’quiet1  plasma  machine,  or  Q-machine,  has  been  used  extensively  in 
experiments  for  examination  of  basic  plasma  phenomena.1  Time  indepen¬ 
dent  equilibria  for  the  single-ended  collisionless  Q-machine  have  been  stud¬ 
ied  analytically  by  a  number  of  authors,  notably  MacIntyre,3  and  Kuhn.3 
The  basic  analytic  approach  is  to  assume  a  steady  state  potential  profile, 
monotonic  in  creasing /decreasing  or  single  extremum,  assume  distribution 
functions  for  the  plasma  source,  usually  half-Maxwellian,  then  model  the 
distribution  functions  as  functions  of  the  potential.  It  is  then  possible,  by 
integration  of  Poisson’s  equation,  to  define  a  ’pseudo-potential’  and  obtain 
equations  that  allow  one  to  solve  for  such  quantities  as  the  plasma  poten¬ 
tial,  the  potential  extremum,  and  the  floating  potential  for  the  collector 
(one  also  gets  values  for  plasma  density).  For  a  review  of  the  technique  the 
reader  is  directed  to  Tidman  and  KralL4 

In  this  paper  we  restrict  our  attention  to  the  negative  biased,  single 
maximum  potential  profile,  since  the  simulation5of  this  case  yielded  results 
at  odds  with  theory.  We  will  describe  the  analysis  for  fixed  bias,  although 
floating  bias  is  easily  handled  with  minor  modification. 

The  one  dimensional  model  for  the  single-ended  Q-machine  is  as  follows. 
The  source  for  the  plasma  is  a  hot  plate,  which  we  will  take  to  be  on  the  left, 
and  distance  is  measured  positive  towards  the  right,  and  at  x  =  L,  is  a  cold 
collector  plate;  both  are  absorbing.  The  ions  are  generated  from  neutral 
gas  according  to  the  Langmuir-Saha  law,  and  hence  have  a  half-Maxwellian 
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■v  \  "  V'  v % 


T^r 


During  Phase  2,  |J(t)|  decreases  on  the  average  but  shows 
some  very  pronounced  oscillations. 

Phase  3  '  "electron  trapping  phase”).  This  phase 

is  characterized  by  the  formation  of  a  trapped  electron  popu¬ 
lation  near  HP  (Figs.  6  and  7) ,  which  leads  to  a  small  negative 
potential  dip  at  the  low-potential  edge  af  the  moving  double 
layer  already  observed  in  Fig.  5.  This  event  is  also  indicated 
by  a  little  dip  in  |<J(t)j  . 

Phase  4  (e^fe  tlb;  "trapped-electron  phase").  Throughout  this 
phase,  the  trapped  electron  population  (which  surrounds  an  ec¬ 
centrically  rotating  electron  hole)  persists  and  remains  local¬ 
ized  in  the  region  USxio.3  (Figs.  8  and  9).  Adjacent  to  its 
right  we  observe  a  region  of  fairly  uniform  and  roughly  half- 
Maxwellian  electrons  ("plasma  region") ,  and  a  region  of  signi¬ 
ficant  electron  acceleration  to  the  right  ("CP  sheath  region"). 
During  all  of  Phase  4,  the  plasma  region  spreads  to  the  right, 
while  the  CP  sheath  recedes. 

The  potential  profile  always  exhibits  a  relative  maximum 
near  HP  and  remains  fairly  flat  in  the  plasma  region,  the  po¬ 
tential  minimum  being  negative  over  most  of  Phase  4  but  going 
positive  at  t»  By  that  time  ,  all  of  the  old  ions  have  ex¬ 
ited  the  diode  to  the  left,  at  a  speed  about  ten  times  their 
thermal  injection  velocity. 

At  t  *  b,  the  cycle  starts  all  over  again. 
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At  t  -  i2  (i.e.,  approximately  at  the  time  the  old  ions  reach 
CP) ,  the  initial  monolayer  near  HP  (Fig.  4)  turns  into  a  dou¬ 
ble  layer  (as  indicated  by  an  inflection  point  in  the  potential 
profile.  Fig.  5)  that  starts  moving  to  the  riqht,  leaving  in 
its  wake  a  region  of  low  potential. 
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Phase  2  (e^6t£e-j;  "electron  ring  phase”).  By  t=e-j ,  the  vor 
tex  observed  during  Phase  1  has  closed  into  what  may  be  loosely 
called  an  "electron  phase  space  ring"  (Fig.  4) .  Phase  2  is  char 
acterized  by  the  "thermalization"  and  exiting  to  the  right  of 
this  ring  (Fig.  5) . 


near  future . 

Close  inspection  of 
the  simulation  data11 
has  led  us  to  tentative¬ 
ly  subdivide  the  full 
cycle  into  four  phases, 
referred  to  as  Phases  1 
through  4  and  described 
below.  Each  of  these 
phases  exhibits  a  char¬ 
acteristic  behavior  of 
electron  phase  space  and 
also  becomes  manifest  in 
the  external  current 
history.  In  Fig.  1,  times 
a,  e-j ,  &2>  e3»  and  15  correspond  to  characteristic  events  in 
electron  phase  space,  while  times  i-j  through  i4  correspond  to 
characteristic  events  in  ion  phase  space.  The  horizontal  line 
marked  Jes  indicates  the  electron  saturation  current,  i.e.,  the 
largest  time-averaged  current  that  could  possibly  be  drawn 
through  the  diode.  For  the  present  run,  Jes  *  -127.7. 

In  what  follows,  each  phase  is  briefly  described  and  illus¬ 
trated  by  two  sets  of  simulation  snapshots  -  one  taken  at  the 
beginning  and  the  other  taken  toward  the  end  of  the  phase  in 
question.  Each  set  contains  a  potential  plot,  an  electron 
phase-space  plot,  and  an  ion  phase-space  plot. 

Phase  1  (a4t4ei;  "electron  vortex  phase").  Conditions 
right  before  the  onset  of  Phase  1  are  almost  identical  with 
those  prevailing  at  the  end  of  Phase  4  {Fig.  9).  By  definition. 
Phase  1  commences  as  the  trapped  electron  population  near  the 
hot  plate  suddenly  separates  from  the  "beam"  electrons,  so  as 
to  form  a  "phase-space  vortex"  (Fig.  2) .  Moving  to  the  right, 
this  vortex  widens  in  both  x  and  v,  and  by  the  end  of  Phase  1 
evolves  into  an  electron  ring. 

During  the  present  phase,  the  space  potential  flips  from 
low  to  high,  the  external  current  magnitude  |J(t)|  goes  through 
a  sharply  peaked  maximum,  and,  at  t=i1 ,  the  initially  coherent 
ion  population  splits  into  the  "old"  ions  (accelerated  to  the 
left)  and  the  incoming  "new"  ions  (immediately  returned  to  HP) . 


tribute  well  beyond  the  present  state  of  understanding  of  the 
PRI. 


2.  MODEL .  The  model  considered  is  the  one-dimensional,  col¬ 
lisionless,  electrostatic  single-emitter  plasma  diode  described, 
e.g.,  in  Refs.  7-9.  Electrons  and  ions  emanate  constantly  from 
the  hot  plate  HP  (x=0,  po-ential  V*0,  temperature  T)  with  half 
Maxwellian  velocity  distributions  corresponding  to  the  temper¬ 
ature  T  and  the  emission  densities  ngQ  and  n^Q,  respectively. 
Particles  hitting  an  electrode  are  absorbed. A  constant  bias  Vc 
is  applied  to  the  cold  plate  CP  (x=L) . 

The  d.c.  equilibria  of  this  diode  are  conveniently  classi¬ 
fied  in  terms  of  the  (»^c,oc)  parameter  diagram,®'9  where  yjc  = 
eVc/kT  and  OC =  n£0/ng0.  A  quantitative  and  detailed  theoretical 
study  of  the  four  most  relevant  parameter  regions  (negative- 
bias  electron-rich,  negative-bias  ion-rich,  positive-bias  elec¬ 
tron  rich,  positive-bias  ion -rich)  was  given  in  Ref.  9,  and  a 
preliminary  comparison  with  pertinent  particle  simulations  was 
presented  in  Ref.  10.  In  the  meantime,  more  systematic  and  de¬ 
tailed  simulation  data  for  all  four  parameter  regions  have  be¬ 
come  available  in  the  form  of  a  16-mm  movie,11  from  which  the 
results  shown  below  have  been  taken.  As  an  example  of  strongly 
unstable  behavior  we  have  picked  out  a  run  whose  nondimension- 
al  parameters  are  located  in  the  positive-bias  ion-rich  region. 
The  simulations  were  all  performed  with  the  new  bounded-plasma 
code  PDW1 . 1 2 

3.  SIMULATION  OF  UNSTABLE  DIODE  BEHAVIOR  IN  THE  POSITIVE- 
BIAS  ION-RICH  REGIME.  In  this  section  we  summarize  the  results 
of  a  simulation  run  characterized  by  the  nondimensional  param¬ 
eters  »jc  =  60,  OC*  4,  A*  50.6  (where  A*  L/  J/^kT/n^e2*)  ,  and 
mi/me  "  40.  More  information  concerning  this  run  can  be  found 
in  Sec.  4a  of  Ref.  10. 

External-current  history.  Fig.  1  shows  the  external-circuit 
current  J(t)  for  the  initial  transient  (fill-up  phase)  and  the 
first  1.5  cycles  of  the  final  dynamic  state,  which  is  strictly 
periodic  with  period  At  *  14.4.  In  what  follows,  we  do  not  con 
aider  the  initial  transient  but  rather  concentrate  on  the  cyc¬ 
lic  behavior.  Moreover,  the  present  account  is  essentially  de¬ 
scriptive;  our  results  are  still  in  the  stage  of  discussion, 
and  a  more  in-depth  analysis  will  hopefully  be  available  in  the 
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Single  Ended  Plasma  Device,  Unstable  States 
PARTICLE  SIMULATIONS  OF  THE  UNSTABLE  STATES  OF  A 
COLLISIONLESS  SINGLE-ENDED  PLASMA  DEVICE 

P.L.  Gray,  S.  Kuhn,**  T.L.  Crystal,  and  C.K.  Birdsall 

Electronics  Research  Laboratory,  University  of  California, 

Berkeley,  CA  94720 

ABSTRACT .  One-dimensional  particle  simulations  of  the  col¬ 
lisionless  single-emitter  plasma  diode  with  a  positively  bi¬ 
ased  collector  plate  are  reported  for  "ion-rich"  operation. 

The  simulations  qualitatively  recover  several  important  fea¬ 
tures  of  the  "potential  relaxation  instability"  observed  ex¬ 
perimentally,  including  "moving  double  layers"  followed  by  a 
current-limiting  potential  dip.  In  addition,  the  simulation 
results  provide  very  detailed  insight  into  the  electron  and  ion 
phase-space  dynamics,  which  can  be  correlated  with  the  exter¬ 
nal-current  history. 

1 .  INTRODUCTION .  It  has  long  been  known  from  both  experi¬ 
ments1”6  and  simulations7  that  a  low-density  single-ended  plas¬ 
ma  device  (such  as  Q-machine  or  thermionic  converter)  may  ex¬ 
hibit  a  large-amplitude,  low-frequency  longitudinal  instabili¬ 
ty  when  the  cold  plate  is  biased  positively  with  respect  to  the 
hot  plate.  For  some  time,  this  instability  was  believed  to  be 
of  the  "standing  ion-acoustic  wave"  type, 3*4  but  more  recent 
experiments  with  improved  potential  diagnostics  have  revealed 
a  much  more  complex  dynamic  behavior,  including  periodic  relax¬ 
ation  from  a  state  of  high  space  potential  to  one  of  low  space 
potential,  via  a  "moving  double  layer"  followed  by  a  current- 
limiting  potential  dip.^  Accordingly,  the  phenomenon  in  ques¬ 
tion  is  now  quite  common) y  referred  to  as  "potential  relaxation 
instability  (PRI)".5'6 

In  the  past,  both  experiments  and  simulations  have  primarily 
concentrated  on  establishing  the  potential  profiles  and/or  car¬ 
rier  density  distributions  associated  with  the  PRI.  To  our  know 
ledge,  however,  no  detailed  investigation  into  the  related  e- 
lectron  and  ion  phase-space  behavior  has  been  reported  to  date, 
and  it  is  mainly  this  aspect  where  we  believe  our  work  to  con- 

♦Paper  presented  at  Second  Symposium  on  Plasma  Double  Layers  and  Related 
Topics,  Innsbruck,  Austria,  July  5-6,  1984. 

♦♦Permanent  address:  Institute  for  Theoretical  Physics,  University  of 
Innsbruck,  Austria. 


observed  experimental¬ 
ly  in  Ref.  2.  The  pe¬ 
riod  is  £ts0.8,  and  the 
fluctuation  level  (de¬ 
fined  as  peak-to-peak 
difference  over  aver¬ 
age  value)  is  about 
1005.  The  evolution 
of  the  potential  pro¬ 
file  during  one  cycle 
of  oscillation  (Fig. 

4)  is  qualitatively 
similar  to  the  one 

a  neaative  potential 


shown  in  Fig.  1b  of  Ref.  2.  In  particular, 
barrier  limiting  the  (electron)  current  builds  up  during  part 
of  the  cycle.  Note  that  the  potential  always  goes  positive 
right  in  front  of  the  hot  plate. 


4b.  Region  1,  parameter  point  (nr  =  l 50  ,a  =  4 )  .  With  the  cold-plate 
bias  this  high,  the  large-amplitude  oscillations  are  no  longer 


pesent,  and  the  diode  exhibits 
an  almost  time-independent  mon- 
otonically  increasing  potential 
as  shown  in  Fig.  5.  This  effect 
of  "high-bias  stabilization" 
has  also  been  observed  experi¬ 
mentally.^  In  the  present  ex¬ 
ample,  there  are  no  ions  for 
and  the  system  does  not 
develop  a  quasi-neutral  plasma 
region.  Still,  the  "plasma" 
properties  predicted  in  Ref.  1 


are  relevant  in  that  they  apply  to  the  inflection  point  of  the 


potential  profile.  E.g.,  the  inf  lection -point  potential  pre- 


th 


dieted  is  =3.37,  in  good  agreement  with  Fig.  5. 


Sa.  Region  2,  parameter  point  (’>=10,3  =  0.1).  Again  we  observe 
coherent,  large-amplitude  oscillations  similar  to  these  of 
Fig.  4  and  Ref.  2.  The  time  evolution  of  the  potential  profile 


PDW1  ,  assuming  a  mas.s  ratio  m./m  =40.  In  each  run  the  diode 
started  out  empty,  went  through  a  “fillup"  phase,  and  eventu¬ 
ally  settled  down  at  its  "final"  state.  The  results  presented 
below  are  all  related  to  the  respective  final  states.  Poten¬ 
tial  profiles  are  displaced  in  the  form"1  n  =  eV/kT  vs .  c,  = 

x//e  kT/n  +  e2  ,  and  time  is  normalised  to  the  ion  transit  time: 
o  eo 

T  =  t/2kT/m.  /  L. 

x 

2.  Region  1  ’  ,  parameter  point  ( n  ~  =  —  3  ,  ct  =  0 . 1  )  .  The  final  state 
is  practically  time-independent  and  exhibits  the  monotonical- 
ly  decreasing  potential  profile  expected  (Fig.  2) ,  the  plasna 
potential  being  in  remarkable  agreement  with  the  theoretically 
predicted  value"5  of  =  -4.40. 


3.  Region  2,  parameter  point  (nc  =  -4  ,ct  =  4 )  .  In  its  final  state 

the  diode  is  again  stable,  exhibiting  a  one-maximum  profile 

as  predicted  (Fig.  3).  However,  the  simulation  profile  has  a 

maximum  potential  ns  =  0.25  and  a  olasma  ooter.tial  r®  =  -0.23 

m 

which  values  differ  significantly  from  the  respective  theoret 

ical  ones,  0.43  and  n"'n=  -0.60.  This  cuar.oitative  disa- 

m  p 

greement  has  been  found  to  be  due  to  the  existence  of  a  trap¬ 
ped  electron  population  ir.  the  region  of  the  potential  maxi¬ 
mum,  which  originates  from  the  more  or  less  turbulent  fillup 
phase  and  is  not  accounted  for  in  Ref.  1. 


4a.  Region  1,  parameter  point  ( "y  =  6 p , a  =  4 )  .  t-  this  case,  the 
final  state  is  characterized  by  coherent  large-amplitude  os¬ 
cillations  similar  to  the  "potential-relaxation  instabilitv" 
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Abstract .  Particle  simulations  of  longitudinal  phenomena  in 
the  collisionless  single-ended  Q  machine  recover  some  bas¬ 
ic  equilibrium  properties  predicted  theoretically,  as  well  as 
the  "potential-relaxation  instability"  and  its  "high-bias  sta¬ 
bilization"  observed  experimentally.  The  results  provide  a 
starting  point  for  more  detailed  investigations  in  various 
directions . 


1.  Introduction .  The  model  considered  is  the  one-dimensional, 
collisionless,  electrostatic  single-emitter  plasma  diode  de¬ 
scribed,  e.g.,  in  Ref.  1.  Electrons  and  ions  emanate  constant¬ 
ly  from  the  hot  plate  (x=0,  potential  V=0,  temperature  T)  with 

half  Maxwellian  velocity  distributions  corresponding  to  the 

_  ....  +  .  + 
temoerature  T  and  the  emission  densities  n  and  n.  ,  resoec- 
*  eo  10 

tively.  Particles  hitting  an  electrode  are  absorbed.  A  con¬ 
stant  bias  Vc  is  applied  to  the  cold  plate  (x=L) . 

In  Ref.  1,  the  d.c.  equilibria  of  this  diode  were  studied 
with  the  assumption  that  there  are  no  trapped  particles  in  the 


potential  wells.  They  were  classified  in  terms  of  the  (r\  ,a) 


parameter  diagram  shown  in  Fig. 

1  ,  where  n  =  eV  /kT  and  ct=n+  /n+  . 

c  c  10  eo 

In  the  present  paper,  we  describe 
simulation  runs  for  four  param¬ 
eter  points  (nc,cO  representa¬ 
tive  of  the  parameter  regions 
1',  2',  1,  and  2.  The  results 
are  compared  with  predictions 

of  Ref.  1  and  recent  exoeriment- 

.  .  t  2,3 

al  ooservaticns . 

The  simulations  were  performed 


_ .  ,  _  ^  .  .  ,  ...  with  the  new  bounded-plasma  code 

Fia . 1  D.c.  potential  profiles  r 


♦Paper  presented  at  International  Conference  on  Plasma  Physics,  Lausanne, 
Switzerland,  June  27-July  3,  1S84. 
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Figure  2.,  the  assumed  potential  profile  for  lingle-maximum  case,  and 
electron  and  ion  distribution  functions  at  various  positions  along  the  pro¬ 
file;  for  electrons  *  *■  =  —  *icp *  =  i^«ui 

*-  ~ V# 3  *  “  tfcp  |  for  lOnS  fi  —  TTl€/M t)  ^io  s  —  i 

•'•a  *  v^:,  v*  =  y/-mp*  ViL.  *  yf-i^Ur 


the  distribution  functions  for  the  two  species  at  different  positions  along 
the  potential  profile,  as  indicated.  The  high  velocity  cutoff  in  f~  merely 
reflects  the  fact  that  electrons  with  initial  energy  greater  than  c<f>ep  make 
H  across  the  system  and  impinge  on  the  collector  plate  (i.e.  the  negative 
half  of  the  distribution  is  due  entirely  to  reflected  electrons).  For  r?  <  0 
(x  >  x3),  we  have  for  the  electron  distribution  function, 

/.(i'.'Kx))  =®cp(-(i/2  —  »|(x)l)  •  B(u  +  yJrj(x)-rtcP)- 

where  again  the  high  negative  velocity  cutoff  indicates  the  electron ‘with 
enough  energy  to  reach  the  cold  plate. 

For  the  ion  distribution  function  we  have,  for  x  <  x\ 

fi{vMx))  =  <*  exp(-(^  +  n(x)))  •  +  y/rj^-ri). 

The  incoming  ions  are  initially  slowed  by  the  potential  bump  and  those  with 
energy  less  than  are  reflected.  For  x  >  xx  we  have, 

UiyMx))  =  a  exp(-(~  +  »7(x)))  •  #(— =  ~  V*lmaz  -  V) 


Those  ions  which  are  not  reflected  by  the  potential  bump  are  accelerated 
towards  the  waff. 


We  would  like  to  obtain  a  set  of  equations  which  will  allow  us  to  cal¬ 
culate  the  potential  maximum  and  the  plasma  potential,  given  a  choice  of 
a  and  ffc,.  The  first  equation  we  will  obtain  from  Poisson’s  equation  and 
the  fact  that  the  electric  field  in  the  plasma  region  is  sero.  We  can  write 
Poisson’s  equation  in  normalized  quantities  as 


dx2 


Qyjf 

=  ".(x^)  -  MX,*?)  s 
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the  quantity  v  is  often  referred  to  as  the  Sagdeev-potential  or  pseudo- 
potential,  and  is  proportional  to  the  square  of  the  electric  field.  Multiplying 
this  by  drj/dx  and  integrating  with  respect  to  x  we  get 


where  C  is  an  arbitary  constant.  We  know  that 


dx 


dx 


—  0 


so  that 


'X«Xl 

•X»X« 

II 

Integrating  the  expressions  for  n,  and  ne  with  respect  to  potential  from 
*7  =  *7m«x  to  q  =  »7P1  we  obtain  the  first  relation 

=  exP(r7p)  (l  +  erfy/rip  -  »fcp) 

-  *xp(*l*uix)  (l  +  crfy/nmas  -  tjcp) 

^  ^(^ep)  ^  ^ *lmoz  *?cp  “  *Icpj 

+  2  ^exp(r;mai)erf 

+  a  exp(-r?p)  (l  -  erf  v/wTT^) 

”  fl  Otp(”fymax)  ^  y/Hmaz  *7p^  = 


X 


where  erf(x)  s  /  exp(-y2)dy.  The  fourth  term  in  the  expression  above 
0 

is  due  to  the  assumed  presence  of  the  hole  in  electron  phase  space. 
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The  second  relation  comes  from  the  quasi-neutrality  condition  in  the 
plasma  region  (n,(»7p)  =  ne(r?p)).  By  integrating  the  distribution  functions 
over  velocity  we  get 


These  two  equations  will  allow  us  to  determine  rjmax  and  r?p  given  initial 

» 

choices  of  a  and  cold  plate  bias  rjep.  For  the  case  of  a  floating  cold  plate  the 
additional  equation  required  is  supplied  by  the  condition  that  the  electron 
and  ion  fluxes  to  the  cold  plate  must  be  equal  at  equilibrium. 

From  our  simulation  study  of  equilibria  in  the  single-ended  Q-machine 
we  have  found  that  the  ’hole’  in  electron  phase  space,  for  single  maximum 
potential  profiles,  is  filled  to  a  MaxweQ-Boltzmann  level.  This  can  be  seen 
from  Figure  3  which  shows  the  potential  profile  for  a  typical  simulation 
for  such  a  case,  as  well  as  the  time  average  distribution  function  for  the 
electrons  windowed  across  the  potential  bump  as  indicated.  We  have  also 
plotted  the  expected  Maxwell-Boltzmann  distribution  on  top  of  the  simu¬ 
lation  result  for  comparison.  Assuming  this  to  be  the  case  simplifies  the 
model  for  the  electron  distribution  function,  which  becomes 

b. 

•  «  i 

M*Mx))  =  exp(-(i/3  -i7(x)))  •  E(v+}/ri(x)-ricp)- 

Having  assumed  that  the  electron  phase  space  hole  is  filled,  the  first  of 
the  two  relations  above  becomes 
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V 


Figure  3.,  the  •ingle-maximum  potential  profile  from  a  limulation  with 
o  *  4.0  and  t fa  *=  floating(ie.  «  -1-84),  and  the  time  averaged  electron 
distribution  function  calculated  at  x  -  Xg  u  indicated  on  the  potential 
profile.  , 
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I’M  ~  =  exp(f?p)  ( 1  +  erf  y/qp  -  ifc p) 

~  cxp(r?maI)  (l  +  erf y/rjmaz  ~  Hep) 

■+■  ^XPijIcp)  ~  Hep  ~  \/Tlp~m  Hep) 

+  a  exp(-i7P)  (l  -  erf  v^ma,  -  f?p) 

“  ®  eXp(“IJ 7 max} ^max  tfp'j  ~  ®* 

it  differs  from  the  realation  for  the  initial  model  only  in  the  abscence  of  the 
term  for  the  phase  space  hole. 

Figures  4  and  S  show  the  comparison  of  simulation  results  with  the 
numerical  calculation  of  the  predicted  values  for  r?p  and  ffmax  for  a  cold 
plate  bias  of  tjep  =  -4.0,  as  functions  of  neutralisation  parameter  a.  As 
can  be  seen  the  calculation  with  the  assumption  of  the  electron  phase  space 
hole  being  filled  to  a  Maxwell- Boltzmann  level  agrees  very  well  with  the 
results  from  our  simulation  study. 

The  fact  that  the  electron  phase  space  hole  fills  to  a  Maxwell-Boltzmann 
level  is  as  yet  unexplained.  We  believe  it  is  due  to  trapping  of  electrons  by 
fluctuation  of  the  potential  profile.  In  reasonably  long  simulations  with 
a  =  4.0  and  a  floating  cold  plate,  we  have  observed  the  trappinng  vortex  to 
expand  into  the  system  while  maintaining  the  potential  extremum  predicted 
by  the  above  analysis.  This  might  lead  one  to  think  that  there  might  exist 
an  extremely  long  period  relaxation  process  as  has  been  observed  in  other 
simulations.5.  We  are  currently  examining  this  posibility. 
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Figure  5.  rjm+x  ▼*.  a  from  numerical  calculations  assuming  the  electron 
phase-space  hole  to  he  filled,  empty,  and  resultes  from  simulations. 
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F.  MULTIFLUID  DERIVATION  OF  THE  ALFVl&N  ION-CYCLOTRON 

LINEAR  DISPERSION  RELATION 


Niels  F.  Otani  (Prof.  C.  K.  Birdsall) 


The  linear  dispersion  relation  appropriate  to  the  Alfven  ion-cyclotron  (AIC)  instability 
is  well-known  for  propagation  of  the  wave  along  the  background  magnetic  field  Bo  in  a 
uniform  plasma  (e.g.,  see  Ref.  1).  Here  we  again  derive  the  AIC  dispersion  relation  but 
use  a  “multifluid”  model  instead  of  the  more  usual  Vlasov  equation  as  the  basis  for  the 
calculation.  The  advantage  of  this  approach  lies  in  its  conceptual  simplicity — one  observes 
the  direct  connection  between  features  of  the  wave  responsible  for  producing  the  various 
linear  currents  and  terms  in  the  dispersion  relation. 

Assume  for  the  moment  an  unperturbed  ion  distribution  of  the  form 

f?(vx,v±,9)  =  f°{vx,v±)p°{9) 

with  the  normalizations 

J  dvMdv±f°(v„v±)  =  1 

r  dBp°{9)  =  1. 

Jo 

The  distribution  function  is  thus  uniform  in  space  and  separable  with  respect  to  the  gy- 
rophase  angle  9.  It  is  not,  however,  time-independent,  since  9  is  not  a  zero-order  constant 
of  the  motion. 

By  considering  the  ion  phase  space  fluid  in  the  Vlasov  description  to  be  composed 
of  an  uncountably  infinite  number  of  cold  fluids,  each  with  its  own  characteristic  set  of 
zero-order  quantities  vx,  v±,  and  9  at  time  f,  we  may  write  the  perturbed  ion  current  in, 
say,  the  y-direction  as 

(•A)»i(*i  t)  =  ej  dvxdv±f°{vx,v±)  J  d9p°(9)(novyl{vx,  z,t)  +  n1(v„ux,ff,z,t)vy)  (1) 

In  this  expression  n0  is  the  unperturbed  ion  density,  vy  =  vj.  cos  9,  and  vvl(vx,  z,t)  and 
dvjdvxfotii(vt,  vx,  9,  z,  f)  are  respectively  the  perturbation  y- velocity  and  perturbation 
density  of  the  cold  fluid  element  with  zero-order  quantities  vx,  v±,  9,  and  z  at  time  t.  A 
mathematically  rigorous  derivation  of  this  infinite  fluid  model  exists  and  will  be  presented 
in  a  future  report. 

The  perturbation  density  may  be  expressed  as: 

Q 

ni{v„v±,9,z,t)  =  -n0~zi{vx,vx,9,z,t) 

Thus 

{Ji)vi  =  ri0e  J  du,dt;i/0(uJ,ux)  J  d9p°(9)(vvl(vx,  z,t)  -  vy-~z\(vx,v±,9,z,t)).  (2) 


Now  consider  the  motion  of  ions  in  a  single,  purely  transverse,  circularly  polarized, 
sinusoidal  wave  of  wavenumber  k  and  frequency  ui  propagating  in  the  direction  of  a  uniform 
magnetic  field  Bo-  The  fields  are  described  by 


B (z,t)  =  B0z  +  8i(-xsin(kz  -  ut)  +y  cos (kz  -  ut)) 
u/B 

E(z,t)  =  -~(xcoa(kz  -  ut )  +  y  ain(A:ar  -  ut)). 
kc 

The  ion  equations  of  motion  in  rectangular  coordinates  are 


dv 

dt 


v  x  B(z,  t) 
c 


) 


(3) 


where  m  is  the  ion  mass. 

Linearized  ion  orbits  may  be  derived  from  Eq.  (3)  given  the  zero-order  orbits 


w*o(0  =  ulo  sin(Ot  +  90) 
vyQ (t)  =  v±0coa(Clt  Oq) 

*o(*)  =  zo  +  vx0t 

where  0  =  uic,  =  eBo/mc.  We  obtain  for  small  wave  perturbations,  (real  part  implied): 


fal”-  n  -  **>  (5> 

U*1  (t)  =  cBi/mc  (VxQ(t)  _  t’t,y0(t))  exp(ikzo(t)  -  iut)  (6) 

*i(0  =  ~  n)3  (Vx°W  ~  ivvo(t))exp(ikzQ(t)  -  iut).  (7) 


Subtituting  these  linearized  fluid  quantities  into  Eq.  (2),  we  obtain  (real  part  implied): 


(. Ji)yt  j  dv,dv±fQ(vz,vx)  j  dffp°(9) 


(u>  —  kv, 

ijj  —  kvx  —  n 


(w  -  kvx  -  n)2 


cos#(cos0  +  isin 


in  9)^j  e' 


kz—iujt 


(8) 


Here  we  note  that  a  time-dependent  unperturbed  distribution  function  will  still  yield  a 
perturbed  current  varying  as  e‘fc,-,wt  provided  /  d9p°(9)  cos2  9  and  /  dOp° (9)  sin  9  cos  9  are 
time-independent.  This  is  true  in  particular  for  the  “four-spokes”  distribution  functions 
used  in  many  of  our  simulations.  We  now  restrict  our  derivation  to  distributions  with 
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I 


I  dffp°(9)  sin 5 cos#  =  0 


which  is  true  of  both  four-spokes  distributions  and  gyrophase-independent  distributions. 
Thus, 


The  electrons  may  be  assumed  to  be  a  cold  fluid,  drifting  in  response  to  the  perturbation 
transverse  electric  field  with  the  E  X  B  velocity: 

4 ff(Je)gl  _  4 irenp  c(Et  x  Bp)y 


_  p  ikz—iuit 


Finally, 


(VxB,),  = 

Ampere’s  Law  then  yields  the  AIC  dispersion  relation: 

0  =  + »x)  (z  Z  a£z  n  +  (w  -C-  fl)« )  (l0) 

The  first  term  is  the  free-space  term  and  the  second  term  is  the  electron  E  x  B  con¬ 
tribution.  From  the  infinite  fluid  derivation,  we  know  the  third  term  comes  from  the 
perturbation  transverse  ion  motion,  while  the  fourth  term  arises  from  the  transverse  ion 
current  generated  by  the  inhomogeneous  first-order  displacements  of  ions  along  Bo-  The 
current  appears  because  the  different  fluids,  with  different  zero-order  transverse  velocities, 
experience  different  differential  displacements,  and  therefore  contribute  with  different  den¬ 
sities  to  the  transverse  current  density  at  a  given  location  z.  It  is  this  last  term  which  is 
responsible  for  the  AIC  instability. 
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POTENTIAL  BARRIER  BETWEEN  HOT  AND  COOL  PLASMAS 

Seigi  Ishiguro,  Inst,  of  Plasma  Physics, 
Nagoya,  Japan  (Prof.  C.  K.  Birdsall) 


A  potential  depression  is  observed  between  two  plasmas  with  different 

electron  temperature  in  Q-machine  experiment  [1].  In  order  to  study  this 

problem,  using  the  PDW1  code,  we  have  performed  one  dimensional  particle 

simulations  with  nonperiodic  boundary  conditions. 

The  conditions  employed  inthe  simulation  are  the  following:  Particles 

are  injected  with  half  Maxwellian  distribution  from  each  sidewall  and 

particles  which  collide  with  the  walls  are  absorbed.  Ion  to  electron  mass 

ratios  are  rru /mg  =  4,  25;  ion  to  electron  temperature  ratio  is  T./T  =  1, 

and  hot  to  cool  electron  temperature  ratio  is  T  ,/ T  *  4.  Hot  and  cool 

en  ec 

electron  thermal  velocities  are  v^  =  1,  vtgc  =  0.5.  The  injected  current 
density  from  the  hot  side  for  electrons  is  Jgh  =  4,  the  system  length  is 
L  =  8,  the  timestep  increment  DT  =  0.02,  and  number  of  time  step  in  a  run  is 
NT  =  6400. 

When  the  system  is  initially  empty,  the  plasma  is  turbulent  in  the 
early  stage,  and  appears  to  seek  equilibrium.  No  interface  is  observed 
between  the  hot  and  cool  side.  See  Figure  1  shows  the  ion  phase  space. 

When  the  system  is  initially  full,  formation  of  ion  holes  is  observed 
as  shown  in  Figure  2,  also  ion  phase  space.  These  ion  holes  move  from  the 
cool  side  to  hot  side  with  -  0.12.  This  velocity  is  smaller  than 

Csh  (=vtehVi7  =  °-2)  and  3reater  than  Csc  (svtecVi;  '  °-I)-  [Note  the 

similarity  between  these  holes  (please  space  vortices)  and  the  real  space 
(x,y)  vortices  between  two  fluids  as  formed  in  the  Kelvin  Helmholtz 
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instability,  at  saturation].  The  potential  of  the  plasma  and  right  side 

wall  (floating)  oscillates  with  plasma  frequency;  with  amplitude 

eVKThot  '  modulated  at  a  low  frequency  (like  uj  /20),  as  seen  in  Figures 

r 

3  and  4.  A  spatial  interface  is  not  observed  as  was  seen  in  the  laboratory 
experiments. 

Some  next  step  may  include: 

(a)  Using  <j> ( x )  time  averaged  (over  at  least  t  ,  oscillation)  in  order 

p  l  asrna 

to  see  whether  the  transient  ion  holes  occur  at  potential  dips; 

(b)  Loading  x,v  space  close  to  the  fg(x,v)  expected  from  equilibrium,  com¬ 
plete  with  ion  hole  and  reflected  electrons  at  t  =  0. 

(c)  Imposing  the  expected  q(x)  for  some  time  0  <  T  <  T<.  (no  Poisson  solver) 
where  T,- is  several  ion  thermal  particle  transit  times,  in  order  to  form 
the  ion  hole.  As  t  >  T<-  the  imposed  i(x)  may  be  slowly  turned  off  and 
the  Poisson  solver  turned  on  (suggested  by  fJ.  Otani ) . 
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ONE-DIMENSIONAL  PARTICLE  SIMULATIONS  OF  THE 
COLLISIONLESS  SINGLE-ENDED  Q  MACHINE1) 

S.Xuhn2) ,  C.K.  Birdsall,  T.L.  Crystal,  P.L.  Gray,  Wm.S.  Lawson 

Plasma  Theory  and  Simulation  Group 
Electronics  Rese-trcn  Laboratory 
University  of  California,  Berkeley,  CA  94720 

Longitudinal  phenomena  in  a  collisionless  single-ended  Q  ma¬ 
chine  are  simulated  with  the  new  bounded-plasma  code  PDW1  /!/. 
Simulations  covering  the  four  main  parameter  domains  dis¬ 
cussed  in  Ref.  /2/  are  presented.  In  the  negative-bias  do¬ 
mains,  the  monotonically  decreasing  and  one-maximum  potential 
distributions  predicted  by  collisionless  d.c.  theory  /2/  are 
recovered,  with  some  quantitative  disagreement  arising  from 
the  filling  of  phase-space  holes  via  numerical  diffusion. 

The  positive-bias  domains  usually  exhibit  the  "potential  re¬ 
laxation  instability" ,  along  with  the  moving  double  layers 
observed  experimentally  /3/.  If,  however,  the  cold-plate  bias 
is  raised  above  some  critical  value,  the  potential  relaxation 
instability  disappears,  giving  way  to  the  one-minimum  or  mon¬ 
otonically  increasing  d.c.  potential  distributions  predicted 
theoretically  /2/.  In  these  cases,  the  electron  sheath  in 
front  of  the  cold  plate  occupies  a  significant  portion  of  the 
interelectrode  gap,  which  is  in  qualitative  agreement  with 
the  experimental  findings  reported  in  Ref.  /4/. 

1)  Work  supoorted  by  DOE  contract  DE— AT03-76ET53604 ,  ONR  con¬ 
tract  N0001 4-77-C-0578 ,  and  Austrian  Research  Funds  con¬ 
tract  S-18/03.  Computations  performed  at  NMFECC,  Livermore. 

2)  Permanent  address:  Institute  for  Theoretical  Physics,  Uni¬ 
versity  of  Innsbruck,  A-6020  Innsbruck,  Austria. 

/I/  Wm.S.  Lawson,  User's  Manual  for  the  PDW1  Code,  Plasma  The¬ 
ory  and  Simulation  Group,  E.R.L.,  U.C.  Berkeley,  CA  94720 
(1  983)  . 

/2/  S.  Kuhn,  Plasma  Phys.  T3,  831  (1931). 

/3/  S.  Iizuka  et  al.,  Phys.  Rev.  Lett.  4_8  ,  1  45  (1  932). 

/4/  G.  Popa  et  al.,  Phys.  Lett.  87A,  175  (1982). 
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PARTICLE  SIMULATIONS  OF  THE  LOW-ALPHA  PIERCE  DIODE 


T.L.  Crystal  and  3.  Kuhn+) 

Plasma  Theory  and  Simulation  Group 
Electronics  Research  Laboratory 
University  of  California,  Berkeley,  CA  94720 

The  evolution  of  small  initial  perturbations  in  the  "classi¬ 
cal"  Pierce, diode  /I/  is  simulated  using  the  new  bounded- 
plasma  code  PDW1  /2/.  The  simulations  cover  the  parameter 
range  0  <  a  <  3tt,  where  a  =  u^L/vq.  In  the  linear  regime,  po¬ 
tential  profiles  of  the  dominant  eigenmode  as  well  as  oscil¬ 
lation  frequencies  of  the  next-to-aominant  eigenmode  are  re¬ 
covered  and  shown  to  be  in  excellent  quantitative  agreement 
with  recent  theoretical  results  /3/.  In  the  linearly  unstable 
cases,  the  nonlinear  final  state  may  be  either  a  new,  non-uni- 
form  d-c.  state,  or  a  state  of  large-amplitude  oscillations. 

In  particular,  for  cc  =  1  .  5tt  the  final  state  is  found  to  depend 
sensitively  on  the  details  of  the  initial  conditions. 

This  work  was  supported  by  DOE  contract  DE-AT03-76ET53064 , 

ONR  contract  N0001 4-77-C-0578 ,  and  Austrian  Research  Funds 
contract  S-18/03.  Computations  performed  at  NMFECC ,  Liver¬ 
more. 

+)  Permanent  address:  Institute  for  Theoretical  Physics,  Uni¬ 
versity  of  Innsbruck,  A-6020  Innsbruck,  Austria. 

/!/  J.R.  Pierce,  J.  Appl.  Phys.  1 5 ,  721  (1944). 

/ 2/  Wm.S.  Lawson,  User's  Manual  for  the  PDW1  code.  Plasma  The¬ 
ory  and  Simulation  Group,  E.R.L.,  U.C.  Berkeley,  CA  94720 
(1983) . 

/3/  S.  Kuhn,  Memorandum  UC3/ERL  M83/61 ,  E.R.L.,  U.C.  Berkeley, 
CA  94720  (1983);  to  appear  in  Phys.  Fluids. 
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The  Classical  Pierce  Diode:  Using  particle  simulations  on 
linear  Sc  oonlinear  behavior  and  final  states 
T.L.  Crystal,  S.  Kuhn*,  and  C.K.  Birdsail. 

EECS/ERL  Cory  Hall,  LMiv  of  Calif  Berkeley  CA  94720 

The  classical  Pierce  diode  is  a  simple  1-d  system  of  two 
shorted  metal  plates,  a  cold  beam  of  electrons  injected 
from  one  side  and  a  neutralizing  background  of  rigid  ions. 
While  the  plasma  medium  is  technically  stable,  the  finite- 
ness  of  the  Pierce  system  allows  stahk  and  unstable  opera¬ 
tion.  It  is  usefully  studied  as  an  archetypical  boijndcd- 
plasma  system,  related  c.g.,  to  Q-inachines,  particle 
accelerators,  thermionic  converters,  and  new  high  power 
microwave  sources;  it  can  also  serve  to  generate  examples 
of  nonlinear,  bifurcating  (arid  oscillating)  equilibria,  virtual 
cathodes  and  double-layers. 

New  particle  simulations  of  the  Pierce  diode  have  success¬ 
fully  recovered  many  novel  Linear  phenomena  including  tin 
dominant  linear  cigenmodos  (seen  in  the  internal  electros¬ 
tatic  fields),  and  the  dominant  and  vt ihdondnant  eigenfrx- 
quencies,  (seen  both  in  the  internal  electrostatics  or,  i  in 
the  external  circuit  current,  </«;*(£))•  These  simulation 
results  conform  very  well  to  detailed  predictions  of  a  new 
linear  analysis1. 

The  final  (nonlinear)  state  recovered  can  shew  critic. »1 
dependence  on  initial  (linear  perturbation.)  conditions,  and 
can  be  made  steady-state  (d.c.)  or  periodic-oseiHc'.or)  by 
simply  changing  trie  initial  conditions  by  a  factor  of  10"-4 
or  less.  A  third  class  of  final  state  is  also  possible  which,  has 
oscillations  which  seem  to  be  nonperiodic.  Such  oscilla¬ 
tions  occur  in  both  the  internal  fields  and  the  external 
current,  Jazt(t)',  of  course,  they  me  di_.  ee  "useful"  or  "ucs- 
tructive”  depending  on  the  ch-etr o.iics  u- vice  bring 
modeled  as  a  Pierce  diode.  The  effects  that  e  sternal  circuit 
parameters  have  on  the  linear  recovered  (slmuiul  ion) 
behavior  of  the  device  are  also  found  to  on  as  predicted1. 

1  S.  Kuhn,  Submitted  to  Phys  fluids  19b  i. 


Work  supported  by  DOE  Contract  Dll- A i’03-76ET530b '.  and 
ONR  Contract  N0001't-77-C-0b7:5;  computations  done  at 
NMFECC,  Livermore. 

•  On  leave  from  the  Univ  of  Innsbruck,  Austria. 
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Axial  Equilibria  of  Collisionlcss  Single-ended  Q-machincs: 
particle  simulations  versus  theory  and  experiment. 

S.  Kuhn*,  P.  Gray.'T.L.  Crystal,  and  C.K.  L)irdsi;!l 
EECS/ERL  Cary  Hall,  Univ  of  Calif  flerkpltty  CA  9-1 79  C 
Only  recently  has  Q-machine  equilibrium  theory  been  sys¬ 
tematically  pushed  to  quantifiable  pri. dictions'.  New  coll;- 
sionless  particle  simulations  have  been  run  with  the  code 
RDW1  in  four  main  operating  regimes  (electron/ ion  r i r - : • 
emission,  and  positive/negative  applied  biases).  For  e.is-j 
where  the  applied  bias  is  negative,  trie  simulations  in  met 
recover  stable  final  states  as  expected,  having  equilibrium 
potential  profile  shapes  that  agree  with  theory.  But  for 
cases  with  ion-rich?emission,  the  simulated  potential  values 
can  differ  from  analytic  predictions  because  equilibrium 
(9/  dt  =0)  theory  does  not  describe  a  system  that  ha« 
evolved ;  collisionless  equilibrium  theory  finds  phase  spam- 
regions  that  are  inaccessible  to  emitted  particles,  yet  simu¬ 
lation  shows  that  these  regions  are  in  fact  filled  in 
Maxwell- Boll rm.ann.ly  during  the  system’s  evolution  {d/ of. 
not  10)  in  reaching  its  eventual  steady  statu.  For  positive 
applied  bias,  the  simulations  generally  cannot  recover  ti.o 
analytically  expected  equilibrium  shapes,  and  in  si  cud  indi¬ 
cate  instability',  whether  the  emission  is  electron  rich  or  ion 
rich;  these  have  character  similar  to  the  low  frequency 
sheath  oscillations  (sometimes  called3  "false  ins»  cbilit.es"> 
and  may  also  relate  to  the  process  of  ’’charge  domain  for¬ 
mation"  in  Gunn  oscillator  theory3. 

Extended  simulations  also  have  suggested  explanations  for 
the  experimentally  observed  (and  overly  mysterious/  h :jh 
positive-bias  stabilization ;  raising  the  applied  bins  high, 
enough  turns  the  simulated  potential  profiles  into  simple 
Child-Langmuir  diode  profiles.  Other  recent  experiments 
have  reported  seeing  a  so-called  moving  double-layer 
(which  might  better  be  re-named);  we  believe  that  these 
may  be  related  to  the  "falsely  unstable"  sheath  oseiliaUen.s3 
already  alluded  to  here. 

1  S.*Kuhn,  Plasma  Phys  ,  23,  1H31  (1931). 

2  N.Rynn,  Phys  rluids,  9,  165  (1966). 

3  M.P.  Shaw,  H.L.  Grubin,  and  P.R.  Soiertion,  The  Cunn- 
Ihlsum  Effect,  Academic  Press  (1979),  p.17. 

Work  ^supported  by  DOE  Contract  D  !•-  A  TO  d  -  7  (>  1 1  !‘b  3  0  G  •*.  and 
ONR  Contract  N0001A-77-C-0576;  eonipul  ation-.  dune  at 

NMFECC,  LLNL  Livermore. 
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AN  ANALYSIS  OF  SINGLE-WAVE  PARTICLE  ORBITS  BASED  ON  CONSTANTS  OF 
MOTION  WITH  APPLICATION  TO  SIMULATIONS  OF  THE  ALFVEN  ION-CYCLOTRON 
INSTABILITY 


Niels  F.  Otani 


Plasma  Theory  a ad  Simulation  Group,  Electronics  Research  Laboratory, 
University  of  California,  Berkeley,  CA  94720 


Three  exact  constants  of  the  motion  have  been  found  for  particles  moving  in  a  single,  circularly- 
polarized  wave  propagating  along  a  uniform  background  magnetic  held.  The  constant  corresponding 
to  the  helical  symmetry  of  the  system  fields,  C,  is  of  particular  interest  because  it  permanently 
restricts  each  particle  to  a  relatively  small  region  of  v^-ti ,  space  which  generally  moves  to  higher 
v,  with  increasing  wave  amplitude.  Confinement  of  particles  to  one  of  these  regions  was  verified 
in  single-wave  computer  simulations  of  the  Alfven-ion-cyclotron  instability.  Existence  of  the  three 
constants  of  motion  reduces  the  problem  to  motion  of  the  particle  on  a  two-dimensional  surface  in  v±- 
v,-\p  space.  When  the  wave  amplitude  is  slowly  varying,  the  motion  on  this  surface  is  characterized 
by  the  existence  and  slow  time-scale  motion  of  fixed  points.  At  least  one  stable  fixed  point  always 
exists  and  its  deeply  trapped  particles  exhibit  an  adiabatic  invariant,  namely,  the  action  J.  It  is 
hypothesized  that  particle  trapping,  fixed  point  motion,  and  the  existence  of  J  may  be  responsible 
for  the  appearence  of  an  edge  in  the  v. -distribution  function  shortly  after  wave  saturation.  The  edge 
in  the  distribution  function  is  also  observed  in  some  multiwave  simulations,  suggesting  that  some  of 
the  single-wave  features  of  the  particle  motion  may  still  be  present  even  though  the  principal  basis 
for  the  analysis,  the  conservation  of  C,  is  no  longer  valid. 


Work  supported  by  the  U.  S.  Department  of  Energy. 
Computations  facilities  provided  by  NMFECC. 
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Beam-Driven  Drift  Instability  in  a  Thermal  Barrier  Cell 
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ABSTRACT 

A  low  frequency  instability  (  -j  «  f>a  )  is  sometimes  seen  in 
the  TMX-U  device  at  1.LN1.  when  the  neutral  beam;  arc  turned  on 
in  the  thermal  barrier  cells.  One  possible  explanation  is  that  a 
drift  wave  couples  with  beam  modus  from  the  injected  ions  in  an 
unstable  fashion.  Cur  previous  work  for  this  instability  (  Bull.  Am. 
Phys.  28  ,  Num.  0.  7Q5  (1901).)  has  bcm'i  lim  *">:  1  •  the  hv-i! 
approximation  for  both  linear  theory  and  partiv'c  -  •  '  lations.  VC. 
have  now  extended  both  the  linear  t  henry  arid  ie  j.  ■.  r!  >hmi!  .- 
tions  to  include  a  more  complete  it mo- 1  • 1  !)•■<•.! 

approximation  is  r>o l  made.  A  compirisoi  -  .  :  :  !••.  -ults 

obtained  from  the  local  approximation  and  the  results  obtained 
from  using  a  nonlocal  model.  Even  wh.cn  the  nonlocal  and  local 
growth  rates  are  very  close  the  instability  rnuy  saturate  m 
different  fashions.  Work  supported  in  part  oy  DOF,  and  in  part  by 
ONR. 
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ERRATA 


Errata  for  Plasma  Physics  via  Computer  Simulation  by  Birdsall  and  Langdon: 

Chop.  2.  p.  10.  In  problem  2-2b,  replace  Purus.  1979  by  Poms.  1978. 

Chap  2.  p.  U.  In  (8),  both  expressions  (r— Ar)/2  should  read  r— A/'  2. 

Chap.  2.  p.  IS  8elow  (12)  ihe  reference  should  be  (Hons,  1970b). 

Chap.  2.  p.  17.  In  the  third  line,  replace  'Chapter  4,  Section  4-14"  with  'Appendix  D.' 

Chap.  2.  p.  79  Problem  2-5c  should  read  and  NG/2—  1  independent  values  of  ImG(A') 
(sine  coefficients)' 

Chap  J.  p.  3S.  In  (3)  change  E  to  £oU1.  : 

Chap.  2.  p.  SO.  At  the  beginning  of  Section  3-14,  delete  'every  time  step.' 

Chap.  4.  p.  59  Above  (7),  change  citation  to  t Bans,  1970b). 

Chap.  4.  p.  65.  At  the  beginning  of  Section  4-6,  'Figure  2-6a.  b.  c'  should  be  'Figure  2-6a.  b.' 

'Figure  2-6b,  d'  should  be  'Figures  2-6a(b),  2-6btb).'_ 

Chap.  4.  p.  63.  In  4-7(2).  the  factor  Z'  should  be  Z'(W(v2Jtv,)). 

Chap.  4.  p.  74.  Replace  E}  everywhere  by  i  £,  1 2.  Below  (2),  read  'As  pCr),  p(x)  ..'  Below 
(3).  the  phrase  "If  we'  should  begin  a  new  sentence.  In  the  seventh  line  from  the  bottom, 
replace  'the  q  P  calculation'  by  '(he  /«A  calculation.' 

Chap.  S.  p.  37.  The  last  line  should  continue  '...and  skip  the  ion  motion,  or  make  m,»me  (or 
<ip«<o p)\  then  correctly  ...' 

Chap.  S.  p.  92.  below  (5)  remove  the  cosier  in  the  expressions  for  .vt  and  v,j. 

Chap.  S,  p.  93.  In  Problem  5-5b,  the  first  reference  should  be  Taylor  ami  McNamara  (197|). 
Chap  5.  p.  94  Change  citation  to  Wuo/f  0949). 

Chap  5,  p.  109.  Figure  5-9c(d).  <M.v)  axis  numbering  should  be  0.3.  0  2.  0.1,  0,  -0.1,  -0.2. 

Chap.  S.  p.  US  Below  (5).  change  citation  to  (Lee  and  Birdsall.  1 979a.  o) 

Chap.  A.  p.  136  In  the  first  paragraph,  change  citation  to  Godfrey  anil  Lanydon  (19'6). 

Chap.  6.  p.  133.  In  ihe  first  paragraph,  change  citation  to  Boris  (1970b) 

Chap  7.  p.  ISO.  In  12).  change  six  to  .v. 

Pan  Two.  p.  IS4  In  the  last  paragraph,  replace  'is  essentially'  with  'is  drawn  ,n  part  from  " 

Chap.  3.  p.  162.  In  (7),  the  summand  should  be  P(k). 

Section  3-S.  p.  162.  3  lines  above  (2).  replace  the  "or"  in  "...  form  or  smoothing  ..."  with  'o|.' 
Chap  3.  p.  173.  Below  (13),  change  8-7(8)  to  8-7(4). 

Chap.  3.  p.  131.  In  Problem  3- 1 3b,  remove  the  squares  from  the  u»'s.  In  Problem  s-i3e. 

change  'this  instability'  to  "the  instability  of  this  Section .' 

Chap.  9.  p.  200  In  the  second  paragraph,  change  the  citation  to  Bum-man.  1967. 

Section  10-1.  p  213.  Add:  Tins  chapter  draws  heavily  on  LuuyJon  119731.’ 

Chap  12.  p.  265.  In  Problem  l2-3d,  /'  - ,  should  be  j'  '■*  i 

Chap  13.  p.  301.  The  citation  Tcireaull.  1933  should  be  Berman  el  at.  1983 

Chap.  14.  p.  343.  Equation  l5)  was  derived  using  (6)  which  does  not  always  hold.  The  poicnual 

due  lo  external  charge  must  be  handled  separately,  see  Section  i0-4  and  Oecyk  (1982). 

Chap  IS.  p  353.  In  the  last  line,  replace  'second  next  section"  with  "Section  15-8'. 

Chap.  IS.  p.  360.  In  the  paragraph  beginning  'Methods  for  ...*,  replace  'discussed  in  ihe  next 
Section*  with  "discussed  in  Section  15-8." 

Chap  IS.  p  366.  In  l3bi.  sin  cklt,  2  m  the  E"  term  should  be  souk  St  In  paragraph  preced¬ 
ing  problem,  replace  nUv)  with  n(k) 

Chap.  IS.  p  3~ l  Above  (10),  change  the  citation  to  Lmdman  1 1 97 5 ) 

Sciuon  15-16.  p.  331.  The  citation  in  die  3rd  paragraph  should  be  Brackhdl  and  Forshm.t.  11985) 
Chap  if.  p.  333.  Above  (a),  ihe  citation  Katamuna  (1980)  should  ho  Katanuma  1 1981 1. 

Chap  16.  p.  407  All  four  o!  the  <r  terms  should  be  instead  of  ' r 

Chap  16.  p.  4/3.  The  parenthetical  comment  at  the  end  of  the  second  paragraph  should  appear 

after  'vanishes'  in  the  previous  sentence 

Appendix  4.  p.  422.  In  two  places,  change  citation  to  Cooley.  Lews  and  IVeich  1 1970). 

References.  In  the  first  printing,  add  the  following  references: 

Oenavit.  J  .  'Pitfalls  in  Particle  Simulations  and  in  Numerical  Solutions  of  the  Vlasov  Equation,' 
Methwlen  and  I'ernaUen  de  Mathemanschcn  Pftvuk.  211,  Peter  Lang,  247-2i9,  1980 
Skflllcrmo.  A  ,  'A  Better  Difference  Scheme  for  the  Laplace  Equation  in  C'..  Imdric.il  Coordi¬ 
nates,'  J.  ('ampin.  I’1'..  47.  160-163,  1932 

Skflllermo,  A.,  and  G  Skftllernio.  'A  Fourier  Analysis  of  Some  Difference  Scr.emcs  'fir  r>c 
Lapiace  Equation  :n  a  System  of  Rotation  Symmetry.'  J  Cmnrut.  Pins.  28,  103-1  14  iyS 
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SECTION  II:  CODE  DEVELOPMENT 


A.  PDW 1 :  Particle  Simulation  of  Non-periodic  Systems 

fhn.  S.  Lawson 


The  development  and  documentation  of  the  PDW1  code  (for  Plasma  Device 
Workshop)  has  now  been  completed.  The  code  is  debugged,  and  fully  docu¬ 
mented.  It  is  available  to  users  of  the  National  Magnetic  Fusion  Energy  Com¬ 
puter  System  on  line,  and  the  user’s  manual  is  avaliable  from  the  Plasma 
Theory  and  Simulation  group  as  memo  UCB/ERL  M84/37.  Documentation 
includes  descriptions  of  important  variables,  subroutines  and  flow  charts,  as 
well  as  examples  of  input  flies  and  the  resulting  output. 

PDW1  is  a  general  purpose  particle  simulation  code  which  is  capable  of 
handling  non-periodic  boundary  conditions  and  a  simple  RLC  external  circuit. 
A  large  variety  of  input  parameters  exist,  allowing  one  to  specify  the  usual 
parameters  of  periodic  simulations,  plus:  thermal,  drift,  and  cut-off  velocities 
for  the  injected  distributions;  DC  and  AC  external  biases  in  the  circuit;  and 
injected  currents.  The  code  works  in  one  spatial  dimension,  and  either  one, 
two,  or  three  velocity  dimensions.  A  magnetic  field  is  allowed  in  the  two  and 
three  velocity  dimension  cases,  and  it  is  allowed  to  be  at  any  angle  in  the  x-z 
plane  in  the  three  velocity  dimension  case. 

PDW1  and  modifications  of  PDW1  are  already  in  much  use  at  our  Plasma 
Theory  and  Simulation  group.  It  has  been  used,  and  is  being  used,  to  simulate 
many  interesting  problems,  such  as  the  Pierce  Diode,  double  layers,  Q- 
machines,  plasma-wall  interactions,  and  oblique  electron  Bernstein  waves. 
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Full  start  case.  Time  evolution 
the  right  wall  potential,  from 
t  ■  80  to  t  *  96. 

Figure  4 

Full  start  case.  Time  evolution 
of  the  plasma  potential,  from 
t  =  80  to  t  a  96. 
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